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Abstract 

We generalize the Hamilton- Jacobi formulation for higher order singular systems 
and obtain the equations of motion as total differential equations. To do this we 
first study the constrains structure present in such systems. 
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1 Introduction 

The study of singular systems has reached a great status in physics since the 
development by Dirac |^, |3| of the generalized Hamiltonian formulation. Since 
then, this formalism has found a wide range of applications in Field Theory 
[0- B H> 01 an d it is still the main tool for the analysis of singular systems. Despite 
the success, it is always interesting to apply different formalisms to the analysis 
of singular systems since they may show new features of the system under study, 
in a similar way to what happens in Classical Dynamics ||. 

Recently an approach based on Hamilton- Jacobi formalism was developed to 
study singular first order systems || This approach consists in using Ca- 
ratheodory's equivalent Lagrangians method to write down the Hamilton- Jacobi 
equation for the system and make use of its singularity to write the equations of 
motion as total differential equations in many variables. This new approach, due 
to its very recent development, has been applied to very few examples QlTJ, [12], [L3|, 
TJj and it is still necessary a better understanding of its features, its advantages 



and disadvantages in the study of singular systems when compared to Dirac's 
Hamiltonian formalism. 

Besides that, theories with higher order Lagrangians (or higher order wave 
equations) are important in the context of many physical problems. Examples 
range from Podolsky's Generalized Electrodynamics |T3| to tachyons (ref. |16| 
and references there in). 

Our aim here is to make a formal generalization of Hamilton- Jacobi formalism 
for singular systems with arbitrarily higher order Lagrangians. This generaliza- 
tion is motivated by the attention that higher order systems have received in 



literature [|T7| , [18|, |19 |. A treatment for the case of second order Lagrangians has 
already been developed by the authors f20fl , but here we will make a more gen- 
eral treatment that will begin with the analysis of the constraints' structure of 
such systems (Sect. 2). Next we will use Caratheodory's equivalent Lagrangians 
method to derive the Hamilton- Jacobi equation for higher order systems (Sect. 
3) and analyze the singular case (Sect. 4). The conclusions will be drawn in Sect. 
5. 



2 Constraints structure in higher order systems 

We will analyze a system described by a Lagrangian dependent up to the K-th 
derivative of the N generalized coordinates g^, i.e. 

L q h qi, q i); q i= 



dt s ' 

were s — 0, 1, K and i — 1, ...N. For such systems the Euler-Lagrange equa- 
tions of motion, obtained through Hamilton's principle of stationary action, will 
be 

K d s I dL 



h dt s 



0) 



0. 



This is a system of N differential equations of 2i^-th order so we need 2KN ini- 
tial conditions to solve it. These conditions are the initial values of gj, g^ ^ 
that describe the velocity phase space (VPS). 

The Hamiltonian formalism for theories with higher order derivatives, that 

(s) 

has been first developed by Ostrogradski pij| , treats the derivatives Q (s = 



0, ...,K — 1) as coordinates. So we will indicate this writing them as Q ;= q( s )i- 
In Ostrogradski's formalism the momenta conjugated respectively to q^K-i)i and 



9(s-i)i ( s — 1? •••) K — 1) are introduced as 



P(K-l)i = —TjcT, UJ 

dL 

P{s-i)i = — P(s)i] s = 1, K - 1. (2) 
Notice that the momenta p( s )j (s > 0) will only be dependent on the derivatives 

(2K-l-s) 

up to 

The Hamiltonian is defined as 

H =2^P(s)i Qi -L\qi,..., Q ij , (3) 

where we use Einstein's summation rule for repeated indexes as will be done 
throughout this paper. Anyway, we will write explicitly the summation over the 
index (s) inside the parenthesis for a question of intelligibility. 
The Hamilton's equations of motion will be written as 

OH 
dH 

P (-)i =-^- = {p (-)ijJ ff} J (5) 
were { , } is the Poisson bracket defined as 

_ K ^ OA OB OB OA 

~ ~ s^o dP(s)i dq(s)i dp(s)i' 

The fundamental Poisson brackets are 

{q( s )i,P(s')j} = 5 ss i5ij; [q( s )i, Q( s ')j} = {p( s )i, P(s')j} = 0, 

were i,j = 1,...,N and s,s' = 0,...,K — 1. With this procedure the phase 
space (PS) is described in terms of the canonical variables and p^ (where 
i — 1, ...,7V and s — 0, ...,K — 1) obeying 27V i\" equations of motion (given by 



equations @ and @) which are first order differential equations. So we have 
to fix 2KN initial conditions to solve these equations. These initial conditions 
are analogue to those needed in Euler-Lagrange equations, but now they are the 
initial values of the canonical variables. 

However, this passage from VPS to PS is only possible if we can solve the 
momenta expressions ([]]) and (0) with respect to the derivatives Q ^ q% ^ so 
that these can be expressed as functions of the canonical variables and eliminated 
from the theory. The necessity of expressing the derivatives 1 $ as functions of 
the canonical variables is clear since these derivatives are present in the Hamilto- 
nian definition (|3|) and must be eliminated from all equations in the Hamiltonian 
formulation. This procedure is completely analogous to the elimination of the 



velocities 9, in the Hamiltonian formalism of a first order system |]5|, 22|. 

The same argument can not be applied to the derivatives q\ , qi ^ since 

(K) 

derivatives higher than q , are not present in the Hamiltonian. But now, the 
necessity of expressing these derivatives as functions of canonical variables comes 
from the fact that they are present in the momenta expressions (0). So, fixing the 
initial conditions of these momenta in the Hamiltonian formulation is equivalent 
to fixing the initial conditions to the derivatives q\ , ''in the Lagrangian 
formulation. The same relation connects the momenta and the derivative 

q i. fixing the initial conditions for the momenta P(k-i)i in the Hamiltonian 
formulation is equivalent to fixing the initial condition for the derivatives q j. 
Then, it is necessary that all the momenta ([!]) and (Q) be linearly independent 

(K) (2K-1) 

functions of the derivatives <? qi so that the latter ones can be solved 

uniquely with respect to the former. 

The expression (|I]) for the momenta P(K-i)i shows that they are dependent 



(K) 

only on derivatives up to Q j, so these derivatives can be solved as functions 

9 i= f(K)i [q( s )f, P(K-i)j) ; s = 0, K — 1, (6) 

if, and only if, the momenta are linearly independent functions of the 

derivatives 1 j. For this, it is necessary that the matrix 

= dp(K-Vi _ d 2 L 

9 q j d q id q j 

be non singular. This matrix H^ is called the Hessian matrix of the system and 
is simply the Jacobian matrix of the change of variables q j— > 

Now, from definitions (|1|) and (|2]), we can see that the momenta P(K-2)i are 
dependent on derivatives up to ^ q\ . In addition, the dependence on ft ' is 
linear and the coefficients are the elements H^ of the Hessian matrix. Considering 
that the derivatives q i have already been eliminated from equations using the 
expression @ the derivatives *" q\ ^ can be solved as functions 

( % )= f(K+i)i [q( s )f p(K-x)j-, v(K~2)j) ; S — 0, K — 1, 

if, and only if, the momenta p(K-2)i are linearly independent functions of the 
derivatives ^ q\ ■ For that it will be necessary that the Jacobian matrix of the 
change of variables ^ q + j — > p(K-2)i, with elements Jy given by 



_ dp {K - 2 )i 

<Ji 



3 .(^+1) 

d qj 
be non singular. 



8 2 L 



{ V -f 9 q\8 q i 

1 u—J(K)u 1 J 



1 u—f(K)v 



1 u—f(K)-u 



Continuing this process, we can use the fact that the momenta p^ s )i are de- 
pendent on derivatives up to ^ qi ^ and, noticing that the highest derivatives 
appear linearly with coefficients that are the elements of the Hessian matrix, show 
that the derivatives ^ can be solved as functions 



(K+p) 

qi -- 



f(K+ P )i (q(s)j; P(K-i)j, ■-, P(K-i-p)j) ; s,p = 0, K — l (8) 



if, and only if, the Jacobian matrix of the change of variables ^ > P(K-i-p)i> 
with elements given by 



J i 



dp(K-i- P )i 



d qj 



("If 



d 2 L 



(K+d) 

1 u=f(K+d)u 



d q id q j 



(K+d) 

1 u=f(K + d)u 



-If H„ 



ij | (K+d) 

1 u=f(K+d)u 



(d = 0, ...,p — 1), is non singular. Consequently, it will be the non singularity of 
the Hessian matrix (|7|) that will determine if the passage from the VPS to PS is 
possible or not. 

Let's suppose now that the Hessian matrix has rank P = N — R. In this 

(K) (2K-1) 

case it will not be possible to express all derivatives qi in the form of 

equation @. Without loss of generality, we can choose the order of coordinates 
in such a way that the Px P sub- matrix in the bottom right corner of the Hessian 
matrix has nonvanishing determinant 



det \\H, 



ab 



det 



3 2 L 



o q a o q b 



^0; a,b = R + l,...,N. 



With this condition, we can only solve the P = N — R derivatives q a as func- 
tions of the coordinates q( 8 )i, the momenta P(K-i)b and the unsolved derivatives 

(K) 

q a (a = 1, R) as follows 

(*) ( (K) \ 

9 a= J(K)a 9(s)ij P{K-l)b] 1 a • 



If we substitute this expression in the momenta definition ([J) for P{K-x)i we 
obtain 



P(K-l)i 



dL 



d q i 



9{K-\)i Q{s)j'i P(K-l)a] 1 a 



(K) 
1 a=f(K)a 



But, since we have P(K-i)a = 9{K-i)ai the other R functions g(K-i) a can n °t 

(K) 

contain the unsolved derivatives Q a or we would be able to solve more of these 
derivatives as functions of the canonical variables, what contradicts the fact that 
the rank of the Hessian matrix is P. So we have the expressions 

P(K-l)a = 9(K~l)a (l(s)j', P(K-l)a) (9) 

which correspond to primary constraints 

${K-l)a = P(K-l)a - 9(K-l)a (l(s)j; P(K-l)a) ~ (10) 

in Dirac's Hamiltonian formalism for singular systems. 

(A"+l) 

Analogously, we can only solve the P derivatives 1 a as functions of the 
coordinates q^i, the momenta P(K-i)b and P{K-2)b and the unsolved derivatives 

(K) (_k- +1 ) 

q a and q a (a = 1, R) as follows 

(K+l) ( (K) (K+l)\ 

q a = J(K+l)a I q(s)U P(K-l)b, P{K-2)b] la, q a I- 

Substituting this expression in momenta definitions @ and using the above 
argument relative to the rank of the Hessian matrix, we obtain for the momenta 
P(K-2)a the following expression 

P(K-2)a = 9(K-2)a \Q(s)j', P{K~l)a, P(K-2)a) 

that corresponds to new primary constraints 

®(K-2)a = P(K-2)a ~ 9{K-2)a {<l{s)j', P(K-l)a, P{K-2)a) ~ 0. 

Continuing this process we find that there will be expressions 

P{K-l-p)a = 9(K~l-p)a{q(s)j', P(K-l)a, P{K-l-p)a) (H) 

(p = 0, K — 1) that correspond to primary constraints 

®(K-l-p)a = P(K-l-p)a ~ 9(K~l-p)a {l(s)j', P(K-l)a, P(K-l-p)a) ~ 0. (12) 



As a result, in a higher order system, the existence of constraints involving a 

given momentum p[K-i)a will imply the existence of constraints involving all P( s )a 

(«) 

momenta conjugated to the derivatives q^ s ) a = Q a (s = 0, K—l) due to the fact 

(K) (2K—1) 

that the derivatives q a ,... } q a can't be expressed as functions of the canonical 
variables. Consequently, if the Hessian matrix has rank P = N — R there are 
KR expressions of the form (|TT| ) that correspond to KR primary constraints in 
Dirac's formalism as given by (|T^|). 

The existence of such constraints ' structure in higher order systems has already 



been noticed by other authors. Nesterenko [17] and Batlle et al. [18] discerned 



the existence of such constraints structure in second order systems, while Saito et 



al. ]19| showed, by different arguments, that the constraints structure exhibited 
above exists for arbitrarily higher order systems. Furthermore, it is important to 
observe that the constraint structure showed above is different for a higher order 
Lagrangian obtained from a lower order one by adding a total time derivative. 
We will not discuss this case here but the reader can find a detailed analysis of 



the constraint structure for such Lagrangians in reference |05 . 



3 Hamilton-Jacobi formalism 

Now we will use Caratheodory's equivalent Lagrangians method to extend the 
Hamilton-Jacobi formalism to a general higher order Lagrangian. The procedure 
described in the sequence can be applied to any higher order Lagrangian and is 
not restricted to a singular one. 



Caratheodory's equivalent Lagrangians method [23] can be easily applied to 



higher order Lagrangians. Given a Lagrangian L qi,Qi, Q i , we can obtain 
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a completely equivalent one given by 

r . T ( (K)\ dS (qi,...,q(K~i)i,t) 

£ ' = H*--"J * ' 

These Lagrangians are equivalent because the action integral given by them 
have simultaneous extremum. So we can choose the function S f<fo, q(K-i)i, 
in such a way that we get an extremum of V and consequently we will get an 
extremum of the Lagrangian L. 

To do this, it is enough to find a set of functions (3( s )i(qj, Q(i)j, 1(s-i)j, t), 
s = 1, K, and S (c[i, q(K-i)i, t) such that 

L' (qi, q(i)i = (3(i)i, • q(K)i = P(K)i, tj = (13) 

and for all neighborhood of g (s) j = /3( s )i(qj,q(i)j, -, 



i ft,-, 9 i >0. 



(14) 



With these conditions satisfied, the Lagrangian L' will have a minimum in 
q<s)i — P(s)i so that the action integral will also have a minimum and the solutions 
of the differential equations given by 

(s) d s qi 



Q(a)i = 1 



dt s 



s = 1, K, will correspond to an extremum of the action integral. 
From the definition of V we have 
T ,_ T ( W\ dS (gj.-.g(AT-i)j.*) ^ q(K-i)j, t) dq (u)i 



u=0 



dq(u)i 



dt 



Using condition fll3|) we obtain 



/ (jf) \ dS (qj, 
L ?j 1 



9* 



K-l 



dS (qj, ...,q {K _i)j,t 



u=0 



9q( u ) 



0, 



Q(s)i—0(s)i 



dS_ 
~dt 



(K) 



K-l 



l U,..., <i , - £ 



U=0 



dS (qj,...,q(K-i)j,t) 
dq( u )i 



9(u+l)i 



Since = (3^ is a minimum point of L' we must have 



9 9j 



1{a)i=P(s)i 



or 



<9L 



dS 



9 ?j 



dt 



Q(a)i—P(a)i 



o, 



dL 



dq(K-i) 

For the same reason we must have 



(K) 

q(s)i=i3(s)i d q i 



1(s)i=P(s)i 



dv 



dL 8 dS 



= 



dL 



d fdS s 



pq{K-\)i dq(K-i)i \dt J 
dS 



Q(s)i=/3( s )i 



_ d 2 S dq {u)j 

_dq {K -i)i dtdq {K -i)i dq(K-2)i ^ dq^ K -i)idq( u )j dt 

dL dS d dS 



3(s)i 



dq(K-i)i dq(K- 2 )i dtdq^K-i)% 



0, 



Q(s)i=P{s)i 



or 



dS 



dq(K-2)i 



1(s)i=fi(s)i 



dL 



d dS 



_dq( K -i)i dtdq( K -i)i_ 



9(s)i—P(s)i 



Following this procedure we have the general expression 



dS 



dq{u-i)i 

were u — 1, K — 1. 



dL d dS 



dq (u)i dtdq {u)i _ 



1{s)i=P(s)i 
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Now, using the definitions for the conjugated momenta given by equations ([!]) 
and (0) in the expressions ( fHf ) and ( |I"8D we obtain 

dS 

P(u)i = a ' u = 0, ...,K - 1. (19) 

OQ(u)i 

So, we can see from equation ( |15|) that, to obtain an extremum of the action, 
we must get a function S uji, q(K-i)i, t) such, that 



(«+i) r / (*) 



where if is 

=51 P(u)i ^"5^ ~ L ( ft; ••! V ?'i ) ( 21 ) 

n=0 V / 

and the momenta are given by equation ([[£]) . 

These are the fundamental equations of the equivalent Lagrangian method, 
and equation (|20|) is the Hamilton- Jacobi partial differential equation (HJPDE). 



4 The singular case 

We consider now the application of the formalism developed in the previous sec- 
tion to a system with a singular higher order Lagrangian. As we showed in Sect. 
2 , when the Hessian matrix has a rank P = N — R the momenta variables will not 
be independent among themselves and we will obtain expressions like equation 
(|Tl~|). We will rewrite these expressions as 

p {u)a = -H (u)a (g (s)i ; p (s)a ) ; u, s = 0, K - 1 (22) 

where we are supposing that the expression for the momentum p( u ) a depends on all 
momenta P( s ) a , although we have showed that the expression for the momentum 
P(u)a is not dependent on any momenta P( s )a with s < u. We do this for simplicity. 
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The Hamiltonian H , given by equation fl2lD , becomes 

n = 2.^ P{n)a Qa +P(K-l)aJ(K)a + 2^ ? a 
u=0 u=0 



— L I q( s )i, Q a , q a — f(K)a \ , (23) 

where a, (3 = 1, R; a = R + 1, iV. On the other hand we have 

dH _ df {K)a dL dL df {K)a _ 

- P(K-l)a— + P(Jf-l)a ^ fx) Jxf ~ U, 

d q a d q a d q a d q a d q a 

so the Hamiltonian Ho does not depend explicitly upon the derivatives q a . 

Now we will adopt the following notation: the time parameter t will be called 
t( s ) = <7( s ) (for any value of s); the coordinates q^ a will be called t( a ) a ; the 
momenta will be called P( a ) a and the momentum pt s \o = -P( s )o will be defined 

as 

P(..o - ^ , (24) 

while i/( s )o = -Ho for any value of s. 

Then, to obtain an extremum of the action integral, we must find a function 
S (t( c )a', q(c)a, t) (c = 0, K — 1) that satisfies the following set of HJPDE 

(dS \ 
t, tf u ) a ] q( u )a'iP{u)a = q = 0, (25) 

(dS \ 
t(u)a',q{u)a',P{u)a = ~R J = 0. (26) 

where s, u — 0, K — 1 and a — 1. R. If we let the index a run from to R 
we can write both equations as 



(d S \ 
t( u )a, q(u)a,P{u)a = q ' J = 0. 

From the above definition above and equation (|23|) we have 

dH{ s)0 _ dL df (K)a K -Hs+1) dH {s) 3 hK)a. n 
OP{u)b q { g> OP{u)b s=0 OP(u)b OP(u)b 



(27) 
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CP(w)6 s =o 0>P(u)b 

where u = 0, if — 1, a = 1, i? and we used the fact that 

% = Q(c+i)i = —jjri c = (J, K - 1. 

Multiplying this equation by (it = dt( s )o we have 

dif( s )0 ^ dH' )a 
d( l{u)b = -jr dt {s)0 + ^ — dq (s)a . 

OP(u)b s=0 OP(u)b 

Using t( s )a = Q(s)a and making the index a run from to R, we have 

K-l Qn< 

dq (u )b=Y, a * Q c %)»- ( 28 ) 
s=0 d P(u)b 

We must call attention to the fact that in the above expression, for a = 0, we 
have the term 

f^b dp {u)b dp (u)b 
that should not be interpreted as 

s =0 d P(u)b OP{u)b 

This somewhat unusual choice of notation allows us to express the results in a 
compact way. 

Noticing that we have the expressions 

K-l Qff' 

dq( u )/3 — —dt( s ) a = 8 su 8 a pdt^ s ) a = dt( s )p 
identically satisfied for a, /3 = 0, 1, i?, we can write the expression (f28|) as 

=E «-^^wa; * = !> ■■■» ( 29 ) 

s=0 °P(u)i 
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If we consider that we have a solution S (qi, q(K-i)h t) °f the set of HJPDE 
given by equation ( p7f) then, differentiating that equation with respect to g( u ) c , 
we obtain 



{s)a , y-^ 0n (s)a 



d 2 S 



#5 



%«) c ^ ^(d)/3 dt( d )pdq( u ) c ^ % d)a dq( d ) a dq(u)c 
for a, /? = 0, 1, R; s, u, d = 0, 1, K — 1 and c = 0, 1, iV. 
From the momenta definitions we can obtain 



(30) 



K-l 



d 2 S K ~ l d 2 S 
dP(u)c = E -JT n & dt (d) p+ -7T ^ dq (d)a . 



(31) 

Now, contracting equation fl30|) with dt/ s ) a and adding the result to equation 
II) we get 



K-l Qfft 

dp(u)c+ E a„ S a dt{ s)a 



K-l 

E 

d=0 
+ 



d 2 S 



K-l 



dq( u )cdq(d)a 
d 2 S 



dq(d) a - E 



(s)a 



s=0 



9p(d)c 



dt 



(s)a 



+ 



If the total differential equation given by (^) are valid, the equation above 



K-l 

E 



dm 



{s)a 



dt 



(s)a 



becomes 



dp( u )c 



^ dm 



E 



(s)a 



dt 



(s)a 



(32) 



s=0 d Q(u)c 

were, as before, u — 0, 1, K — 1; c = 0, 1, iV and a = 0, 1, Z?. 

Making Z = S (t( s ) a ', <?0)a) and using the momenta definitions together with 
equation (^5|) we have 

K-l QCj K-l Qg 

dZ=Y. & dt {d)p + -7T dq {d)a , 

d=0 dt (d)P d=0 d( Kd)a 



K-l K-l /K-l 

dZ '= - E H (d)pdt {d) p+ P(d)a E 

d=0 d=0 \ s=0 



/k ~ i dm 



0)q 



dp( d ) a 



dt 



(s)a 



16 



With a little change of indexes we get 



K-l / K-l QJJ> \ 

dZ = J2[ -H m + J2 H^J^T 1 dt W- (33) 



This equation together with equations ( p29|) and fl32| ) are the total differential 
equations for the characteristics curves of the HJPDE given by equation fl2"7|) and, 
if they form a completely integrable set, their simultaneous solutions determine 
S (t(s)a', Q(s)a) uniquely from the initial conditions. Besides that, equations (|29|) 



and (|32|) are the equations of motion of the system written as total differential 



equations. 



5 Conclusions 

We have obtained the equations of motion for the canonical variables of a singular 
higher order system as total differential equations. Each coordinate q( s ) a = 
(a = l,...,R) is treated as a parameter that describes the system evolution. 
The Hamiltonians H'^ a will be the generators of the canonical transformations 
parametrized by tr s \ a in the same way the Hamiltonian Hq is the generator of 
time evolution. If we have K = 1 the results obtained here will reduce to the 
case of first order systems showed in ref. 0. For K = 2 we have the same re- 
sults obtained for a second order system of ref. p0[] , where the Hamilton- Jacobi 
formalism was applied to Podolsky generalized electrodynamics and the results 
were compared to Dirac's Hamiltonian formalism. 

The integrability conditions that have to be satisfied by equations (P2"|) 
and (|33D are analogous to those that have to be satisfied in the first order case. 
These conditions have been derived in ref.[|Kj and can be easily applied to the 
higher order case developed here. These integrability conditions are equivalent 
to the consistence conditions in Dirac's formalism. 
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We must point out that one of the reasons to consider the constraints' struc- 
ture described in Sect. 2 is the fact that if we had considered only the constraints 
containing the momenta P(K-i)a, given by equations of the form (g), when de- 
veloping the singular case in Sect. 4, the coordinate t^ K _i) a = q(K-i) a would be 
an arbitrary parameter in the formalism but the coordinate t(K-2)a = Q(K-2)a 
(that obeys q^K-i)a =Q{K-2)a) would have a dynamics of its own. So, when we 
choose to deal with all constraints given by expression ( |TTD we are avoiding such 
contradictions. Furthermore, if we had not made this choice, we would have 
an unnecessary extra work when analyzing the integrability conditions since the 
constraints involving momenta P( s )a with s < K — 1 would appear in this stage 
imposing extra integrability conditions. 

As we mentioned in Introduction, Hamilton- Jacobi formalism is not well stud- 
ied for singular systems yet. We still lack a complete analysis of the relation 
between the procedures in this new formalism for singular systems and tradi- 
tional ones, specially the relation with Dirac's Hamiltonian formalism. Besides, 
Hamilton- Jacobi formalism shall be applied to various physical systems so that 
we can get a better understanding of its potential to deal with specific problems. 
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